Quantum Wells in Polar-Nonpolar Oxide Heterojunction Systems 
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We address the electronic structure of quantum wells in polar-nonpolar oxide heterojunction 
systems focusing on the case of non-polar BaV03 wells surrounded by polar LaTiOs barriers. Our 
discussion is based on a density functional description using the local spin density approximation 
with local correlation corrections (LSDA+U). We conclude that a variety of quite different two- 
dimensional electron systems can occur at interfaces between insulating materials depending on 
band line-ups and on the geometrical arrangement of polarity discontinuities. 
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I. INTRODUCTION 



Complex oxide materials exhibit a wide vari- 
ety of exotic phenomena from high temperature 
superconductivity^, to colossal magnetoresistance 2 -, to 
electronic and orbital order—. Recently, inspired by ad- 
vances in epitaxial growth techniques, there has been 
increasing interesls^SiLSiSiiS in making progress toward 
the goal of controlling oxide properties by designing 
layered structures at the atomic level. The class of 
artificial materials which has been studied most ex- 
tensively is ABO3 perovskites. Already there is evi- 
dence for the emergence of superconductivity from non- 
superconducting materials^ and magnetism from non- 
magnetic materials^. One element which distinquishcs 
oxide heterojunction systems from weakly correlated 
semiconductor heterojunctions is the possibility of ex- 
ploiting gradients in polarity! 1 - to help control electronic 
properties. In this article we study the electronic proper- 
ties of oxide quantum wells in which the barrier and well 
materials have very different polarities, using the case 
of non-polar BaV03 wells surrounded by polar LaTi03 
barriers as an example. 

Our paper is organized as follows. In Section II, we 
present a qualitative discussion of four possible polar- 
nonpolar quantum well geometries. We have tested con- 
jectures made in Section II by performing local spin- 
density-approximation with local-correlation corrections 
(LSDA+U) calculations for BVO quantum wells in LTO 
and LTO quantum wells in BVO. Throughout this paper, 
we use LTO as shorthand for the perovskites associated 
with LaO and Ti02 layers and BVO as shorthand for the 
perovskites associated with BaO and VO2 layers. In Sec- 
tion III, we briefly detail the LSDA+U electronic struc- 
ture calculation methods we have employed to model 
LTO/BVO/LTO or BVO/LTO/BVO quantum well sys- 
tems. The combination (LTO/BVO) was chosen because 
LTO is polar and BVO is not and because this material 
combination should not lead to large strains at the inter- 
face. Partly for convenience and partly because multi- 
quantum-well systems are likely to be of equal experi- 
mental interest, we assume periodically repeated quan- 



tum wells so that we have a superlattice in the growth 
direction. In bulk both LTO and BVO are d 1 Mott in- 
sulators. We focus on heterojunction systems composed 
of materials which have the same d-band occupancies in 
the bulk because more possibilities exist for interface elec- 
tronic structure compared to the case of heterojunctions 
between d 1 and d° materials such as LaTiOa/SrTiOs, for 
which interface metallicity appears to be inevitable. In 
Section IV we present the results of our ab initio elec- 
tronic structure calculations for the two representative 
quantum well systems and compare with the conjectures 
we made in section II. Finally in Section V we briefly 
summarize our results. 



II. POLAR-NONPOLAR QUANTUM WELL 
STRUCTURES 



In this section, we discuss four different possible ge- 
ometries for a polar-nonpolar quantum well system. The 
main points we wish to make are shown in Fig. [1] which 
illustrates approximate charge transfers and electric po- 
tentials. Throughout this discussion we use LTO to rep- 
resent a prototypical d 1 polar perovskite and BVO to rep- 
resent a prototypical d 1 nonpolar perovskite. We identify 
two possible geometries for both LTO/BVO/LTO (case 
(a) and case (c))and BVO/LTO/BVO (case (b) and case 
(d))quantum well systems. At an interface between two 
strongly correlated bulk materials we should expect a 
complex interplay between the ionic potential energy, ki- 
netic energy, and correlation energies of the two materi- 
als. By using density functional theory with local correla- 
tions (LSDA+U) we can take into account those energies 
approximately with a single electron picture. The elec- 
tronic structure is characterized most crudely by an ionic 
electron counting picture which neglects hybridization 
between oxygen and transition metal orbitals and ignores 
the finite thickness of the various atomic layers. Using 
this picture, the La +3 0~ 2 planes in LTO are positively 
charged and the Ti^O^ 4 are negatively charged. In 
bulk BVO on the other hand, the Ba+ 2 0~ 2 and V +4 0^ 4 
planes are neutral in the bulk. The transition-metal ions 
in both materials then have one electron in the d orbitals 
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to maintain charge neutrality. When we combine po- 
lar perovskites with nonpolar ones, the transition-metal 
layer will distribute the d-charges in a different way to 
compensate for the nonuniform distribution of ionized 
layers in this system, as illustrated in Fig. [T] The distri- 
butions labeled ideal conjecture in this figure assign the 
charge deficiency or excess to the two transition metal 
layers adjacent to the quantum well whereas the dis- 
tributions labeled realistic conjecture assume additional 
charge redistribution Sp due to the difference between V 
and Ti d-orbital energies. We discuss only cases in which 
the center of the quantum well is a center of inversion 
symmetry. In each case the number of layers which would 
have a positive charge in a bulk polar crystal differs from 
the number of layers which would have a negative charge. 

In case (a) , we have two BaO layers and one VO2 layer 
inserted in the unit cell of a superlattice with six LaO + 
layers but seven Ti02 layers. The charges in the super- 
cell cannot be neutral if we assume that each TiC>2 is 
negatively charged as in the bulk LTO states. Instead, 
we expect that the extra electrons will be distributed 
evenly among the TiC>2 layers near the central VO2 layer. 
The easiest way to accomplish this is to have the TiC>2 
layers next to BaO layers to be Ti +35 02~ 4 rather than 
Ti +3 02~ 4 , i.e. that they contain 0.5 d-electrons per Ti 
rather than 1.0 ci-electrons per titanium. This idealized 
conjecture ignores d—p hybridization as well as the differ- 
ence between the d-orbital energies on Ti and V sites. In 
particular a V +4 ion should have a lower intrinsic atomic 
levels than the Ti +3 ion. We should expect that charges 
might react locally to this difference by transferring ad- 
ditional electrons 5p to the VO2 well from nearby TiC>2 
layers. This additional charge transfer partially screens 
the original atomic energy level difference. These consid- 
erations are illustrated schematically in Fig. [1] The d- 
electron system electronic structure reconstructs in this 
way near the interface in reaction to the electrostatic po- 
tentials produced by the bulk crystal layer polarizations. 
Qualitatively similar considerations would apply to wider 
BVO quantum wells. 

In case (b), we have two LaO + layers and one TiC>2 
barrier layer in the middle of the supercell. By the same 
arguments as case (a), the simplest conjecture is that 
the VO2 layers nearest the TiC>2 layers will be negatively 
charged (V+^O^ 4 ) to keep the system neutral. The 
charge configuration Ti +3 0^~ 4 is maintained by the pro- 
tection of two surrounding LaO + layers. Therefore, there 
will be 1.5 d-electrons on V sites in the V +3 ' 5 C>2~ 4 layers 
and 9 d-electrons in the supercell. However, the local 
physics is different in this case. The electron can respond 
locally by transferring electrons from the barrier TiC>2 to 
the VO2 well layers. The reconstructed potential should 
therefore be lower in the well region in this case. 

In case (c) , we switch the ordering of layers from case 
(a) (BaO <-> V0 2 and LaO <-> Ti0 2 ). There are 7 pos- 
itive charged LaO + layers and 6 Ti^O^ 4 layers. To 
compensate for the charge difference, VO2 layers have to 
be negatively charged as V +3 ' 5 0^~ 4 in the idealized con- 
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FIG. 1: (color online) Approximate charge density, electric 
field, and electric potential distributions in different polar- 
nonpolar quantum well systems with inversion symmetry in 
the quantum well center. The arrows (black on line) indicate 
the expected charge densities p centered on different atomic 
layers, the piece-wise constant thick lines lines (blue on line) 
indicate the expected electric fields E between the atomic 
layers, and the cusped lines (red on line) indicate the expected 
electric potential —V(z) distribution. The signs + and — 
stand for the direction of electric field E, the value for p, and 
the potential —V(z). 



jecture illustrated. Additional charges can transfer from 
adjacent Ti^O^ 4 to the U +35 07 4 layers. The recon- 
structed potential should therefore be higher in the well 
region. 

Similarly in case (d), the structure has been changed 
from case (b) by switching the ordering BaO <-» VO2 and 
LaO «-» TiC>2. There is only one LaO + layer in this case. 
The abutting TiC>2 layers should have the configuration 
rpj+3.5Q-4 ^ Q cnar g e neutrality, resulting in a lower 
reconstructed local potential in the center of the cell. No- 
tice that the magnitude of the additional charge trans- 
fer and its spatial distribution expected in the realistic 
pictures needs to be determined by self-consistent micro- 
scopic calculations. In the next section, we will focus on 
cases (a) and (b) to check the degree to which the qualita- 
tive pictures explained here agrees with the results of mi- 
croscopic DFT calculations and to obtain further insight 
into the nature of the two-dimensional electron systems 
which these electronic reconstructions enforce. To avoid 
confusion to the readers, in microscopic calculation, the 
d electrons should be interpreted as the electrons occupy- 
ing the bands, which hybridized with s, p orbitals, with 
dominant d character at transition-metal sites. 
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FIG. 2: (color online) Half of the inversion symmetric unit 
cell. The quantum well consists of two layers of BaO and 
a monolayer of VO2 and is embedded in bulk LaTiC>3. The 
light-orange shadowed region indicating the narrow active re- 
gion in which the electronic interface reconstruction (EIR) 
driven by the polar discontinuity takes place. Oxygen atoms 
are located at the intersections of the checkered lines and form 
octahedra around the Ti and V atoms. The numbers specify 
labels for different layers. The downward arrow indicates the 
location of the VO2 inversion plane. 



IV. LSDA+U QUANTUM- WELL ELECTRONIC 
STRUCTURE 

Case A : (LaTi0 3 )3/(BaO) 2 V0 2 /(LaTi0 3 )3Ti02 



III. ELECTRONIC STRUCTURE 
CALCULATIONS 

In our ah initio electronic structure calculations we 
have considered symmetric quantum well geometries 
only, use the experimental atom position er 12 ! 13 for LTO, 
and neglect atomic relaxation in the BVO layers. (The 
room temperature crystal structure of BVO is not known 
experimentally.) The electronic structure calculations 
were performed using DFT with LSDA+U as imple- 
mented in the software package VAS P 14 i 15 ' 16 , Projector- 
augmented wave pseudopotentialsii are used to describe 
the electron-ion interaction. We sample the full super- 
lattice Brillouin zone (BZ) with a 5 x 5 x 3 mesh and 
used an energy cutoff of 410 eV which we found to be 
sufficient to reproduce bulk LTO properties. We used 
PAW pseudopotentials supplied by the VASP code for 
La, O, Ba, Ti and V but for Ba, Ti and V, we used the 
one with semicore states treated as valance states. The 
screened local correlation U in the partially filled d bands 
of the transition metal elements were treated using the 
rotationally invariant LSDA+U method due to Dudarev 
et al* 1 ^. Since the La / bands lie much higher in energy 
in experiments^^ we impose a much larger C/ iQ = 11 eV 
on La / states to prevent their mixing with d bands at 
low energy. We chose Uxi = 3.2 eV and Jxi — 0.9 eV 
since these values reproduces the Mott-insulating nature 
of bulk LTO 1 ^ and we take Uy and Jy to be the same 
as U~Ti and Jn- We have varied Uy to check the robust- 
ness of our calculations and find that our conclusions are 
independent of this LSDA+U parameter. 



The unit cell of the superlattice consists of a single VO2 
layer at the center of the well surrounded by two BaO 
layers as illustrated in Fig. [51 Because of the smaller p-d 
energy separation in bulk BVO compared to LTO, this 
material combination should form a quantum well cen- 
tered on the VO2 layer, assuming that the p-bands of 
the two materials are nearly aligned at the interface. In 
our superlattice the BVO quantum well layers are sep- 
arated by 7 Ti02 transition metal layers, each in turn 
separated by a LaO layer. Half of the inversion sym- 
metric unit cell is shown in Fig. [5] in which the VO2 
inversion plane is marked by an arrow. Ti and V in bulk 
LTO and BVO respectively both have nominal valence 
charge d . In the quantum well structure we study here 
the number of LaO layers, which are positively charged 
in the polar LTO crystal, is one smaller than the num- 
ber of Ti02 layers. It follows that the total d charge 
per transition metal summed over all layers is expected 
to be reduced by 1. Unless the inversion symmetry is 
broken we should expect the missing d-charge to be sym- 
metrically distributed around the quantum well center. 
The simplified assumptions of the previous section would 
divide this valence change between the two Ti02 layers 
adjacent to the quantum well. 

Fig. [3] shows the (i-projected DOS at a Ti site and a V 
site in each Ti02 and VO2 plane. The states below —1.5 
eV are the d — p bonding states between Ti, V atoms 
and their neighboring oxygens which have dominant p- 
character, whereas the states near the Fermi level (E = 
Ep) are the anti-bonding states which have dominant 
d-character. As we move through the superlattice from 
TO2 barrier layers to the VO2 we see the DOS evolve 
from the typical bulk LTO form (layers 1 to 3) to a 
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FIG. 3: (color online) Layer-resolved d-projected density of 
states at Ti and V sites versus energy. The layer labeling is 
denned in Fig. [2] 



new form near the heterojunction. In particular, there 
is a large depletion away from atomic-like d 1 -charges in 
layer 4. The electronic reconstruction which accommo- 
dates the band offsets and the polarity discontinuity is 
evident very strongly in the layers 5 and 6 surrounding 
the quantum- well layer 4 (layer 6 is not shown in Fig. [3] 
to avoid redundancies because of the inversion symmetry 
of the system). We also observe that the p-d separation is 
maintained at the VO2 layer as expected after the elec- 
tronic reconstruction. The total electron transfer from 
layer 4 is larger than in the ideal conjecture, as expected, 
and more remote layers participate in the transfer be- 
cause of oxygen mediated coherence between transition 
metal orbitals centered on different layers. 

Most Mott insulators are orbitally and magnetically 
ordered. In order to gain some insight into the way 
in which the magnetic order is disturbed by the quan- 
tum well structure, we plot in Fig. U the spin-resolved 
d-projected DOS in all layers. In our calculations we find 
G-type antiferromagnetic (AFM) order in which neigh- 
boring transition metals ions have opposite spins both 
within and between layers for the bulk like TiC>2 layers, 
(layers 1, 2, 3 in Fig.[3]). Layer 4 on the other hand has an 
overall magnetic moment of around 0.1 /j.b per Ti ion in- 
dicating a ferromagnetic layer. Because of the electronic 
interface reconstruction (EIR), the magnetic order in the 
VO2 plane changes from AFM to ferromagnetic. The 
projected moments on each Ti ions in the AFM bulk-like 
TiC>2 layers is calculated to be 0.74 ijlb inside the Ti ionic 
radius which is fixed in our calculation by the pseudopo- 
tential construction. We believe that by enlarging the 
Ti ionic radius, we can get d-orbital integrated moments 
close to the 1 \ib value expected for d 1 Mott insulators. 
Note that ferromagnetism is induced in the VO2 layer by 
the two adjacent strongly depopulated and ferromagnetic 
TiC>2 layers. 

Fig.[5]shows the superlattice band structure along high 
symmetry directions in the Brillouin zone. The bands 



(five majority (red) spin bands and five minority (blue) 
spin bands), which lie far below the Ej? (in the energy 
range between between —1 eV and —0.6 eV), can be iden- 
tified as bulk LTO d-charge character bands by compar- 
ing with Fig. [4] and calculating band-decomposed local 
charge density contributions. In the supercell, we have 
5 bulk TiC>2 layers (layer 1, 2, 3, 7, 8) in our super- 
lattice corresponding in the antiferromagnetic state to 
10 full bands (5 spin-up and 5 spin-down). Therefore, 
there should be ten full bulk like LTO d-bands, iden- 
tified by the solid circle, in Fig. [5] Two full majority 
spin bands and one full minority spin band (near -1 eV) 
highlighted by two dashed circles are associated with the 
d-charges in the VO2 layers. The two majority bands 
close to the interface states in energy and are hybridized 
along the layer-growth direction and delocalized (This 
interpretation has been verified by separating the charge 
density contribution of those bands) . Because of this 
hybridization, the total d-charges (2.4053 e~) at the V 
sites is smaller than the 3 e~ value which would be im- 
plied by a local interpretation of the band structure. It 
can in general be misleading to estimate the number of 
charges at transition metal ions by counting the number 
of reconstructed bands with a particular dominant orbital 
character. The three partially full bands (one minority 
band and two majority bands) near the fermi surface Ep 
are the bands responsible for the charge depletion at the 
Ti02 layers adjacent to the quantum well (layer 4 and 6). 
The momentum k z dependence of the electronic states is 
shown in the right panel of Fig. [5] (from momentum (0 
, 0, 0) at T to (0, 0, k z ) at Z). It shows that, the states 
at the Fermi energy form a two-dimensional electron gas 
centered on three layers, the active region. Notice that 
there also appears a flat band between M and X, which 
corresponds to the direction away from the bonding di- 
rection between Ti and O atoms in charge depleting the 
Ti02 layers, due to the fact there is negligible tunneling 
of electrons along this direction. The two dimensional 
system is metallic in this case and has three partially 
filled bands. 

Since we have discussed the local potential for the elec- 
tron in case (a) in Fig. [T] based on qualitative model 
in which we assume the charges at each layer are con- 
centrated in a single layer along the layer direction, we 
would like to examine how the the crude model predic- 
tion compares with actual DFT predictions. In Fig. [51 the 
DFT plane-averaged Hartree potential for the electrons 
is shown. First Ti02 layer starts at z = and layers re- 
peat with a spatial period of 3.96 A and with the central 
layer (z = 15.84 A) of the supercell substituted by the 
VO2 layer. First LaO layer is located at z = 1.98 A and 
these layers repeat with a period of 3.96 A except that 
the two layers surrounding the central VO2 layer corre- 
spond instead to BaO layers. By inspecting the variation 
of the local minima in the potential, we can see that the 
qualitative behavior anticipated in case (a) in Fig. [T] is 
valid, as shown by the black-dashed line. The electric 
field in the qualitative model should be considered as a 
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FIG. 4: (color online) Spin and layer resolved densities-of- 
states versus energy. The labels on top are layer labels as 
in Fig. [2] Since we allow G-type AFM order in the barrier 
material, there are two distinct transition metal atoms in each 
layer. The DOS for majority spins and minority spins are 
represented by red and blue lines respectively and the Fermi 
level is marked with a horizontal dark line at E = Ef- 



macroscopically averaged electric field between potential 
minima in the microscopic calculations. Notice that the 
potential of the VO2 layer is somewhat lower than antic- 
ipated. 

In Fig. we show the linear electron number density 
\{z) = 0.5 J p(x, y, z)dxdy along the layer-growth direc- 
tion. The charge number spatial distribution p(x, y, z) 
is normalized per transverse unit cell which is doubled 
because of the system's antiferromagnetic order, and in- 
cludes the total number of electrons contributed by the 
(i-bands between with energies up to 1.5 eV below the 
Fermi level. (See Fig. [5]) The factor 0.5 is used to convert 
the normalization to d-electrons per transition metal ion, 
to facilitate comparison with the qualitative picture. The 
dashed lines represent the boundaries for the transition- 
metal oxide layers used in this electron number construc- 
tion. The area under the curve A(z) represent the total 
number of electrons from all (i-bands. We see an ap- 
parent charge depletion at the TiC>2 layers centered at 
z = 11.88 A and z = 19.8 A along with large charge ac- 
cumulation at the VO2 layer. There is only minor charge 
depletion at second-nearest TiC>2 layers [z — 7.92 A and 
z = 23.76 A) away from VO2 layer. The total charges 
from the (i-bands, 7 e~, is obtained by adding the number 
of electrons at each layers and agrees with the qualita- 
tive predictions for case (a). The number of electrons in 
the depleted TiC>2 layers is 0.3793. The other layers have 
small deviation from one electron. The calculation shows 
a clear deviation from ideal conjecture we made for this 
case and is closer to the realistic conjecture. In general, 
the charges in the well are mainly transferred from the 
adjacent TiC>2 layers but there is also a minor contribu- 
tion from other TiC>2 layers as shown from the number 
of electrons in each layers. 
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FIG. 5: (color online) Spin-resolved band structure at the 
Brillouin zone. The horizontal axis represents superlat- 
tice crystal momentum. The red solid lines represents the 
majority-spin bands and the blue ones are the minority-spin 
bands. The momentum dependence along the layer direction 
k z is indicated in the region between the F and Z points. 
The full circle marks the d 1 bands for Ti atoms in LTO bulk 
states. The dashed circles represents the d bands for V atoms 
in VO2. The bands associated with EIR are flat along k z 
demonstrating their quasi- two-dimensional character. The 
shadowed blue square represents the projection of the super- 
lattice Brillouin zone in the plane of the quantum well. 
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FIG. 6: (color online) Plane averaged Hartree potential along 
the layer-growth direction z. The blue curve represents the 
Hartree potential experienced by electrons — eV(z). The 
dashed line indicates the net potential variation between 
charge-concentrations centers in each layers. 

Case B : (BaV0 3 )3/(LaO)2Ti0 2 /(BaV0 3 )3V02 

In this section, we consider the quantum well structure 
shown in case (b). In this case, the number of LaO layers, 
which are positively charged in the polar LTO crystal, is 
greater than the number of TiC>2 layer by one. It follows 



6 




LayeM /> *J\I\ 


l/\ ^ 


■ Vi TTTTT MM fj. 

Layer 2 A 1 


V 

1 ri l-Nl — r— 


Layers 




— I — I — I — I — h-T^TT 1 T 1 1 I ^jtfH- l — l — 
Layer 4 1 1 


H — 1 — i|H — l> 1 


Layer 5 /I 


H — 1 — 1 — 1 — 1 — 1 — h-T" 
Ti 



3.96 7.92 11.88 15.84 19.8 23.76 27.72 31.67 



FIG. 7: (color online) Linear electron number density X(z) 
along layer-growth direction z. The dashed lines represent 
the boundaries used for the partitioning of density into con- 
tributions from different layers. The total number of electrons 
in each layer contributed by the d-character bands in Fig. [5] 
is indicated by the numerical values between dashed lines. 
In this case, the total number of electrons in the supercell is 
seven. 



that the total d-electrons per transition-metal summed 
over all layers is expected to be increased by 1. Unless 
the inversion symmetry is broken we should expect the 
extra d-charge to be symmetrically distributed around 
the quantum well center. 

Fig. [8] shows the d-projected DOS at a Ti site and a V 
site in each TO2 and VO2 plane. The states below —1 
eV are the p — ci-bonding states between the transition- 
metal atoms (Ti, V) and the p orbitals of the neighbor- 
ing oxygens and have dominant p character. The states 
near the Fermi level (E = Ejr), on the other hand, have 
dominant d-character. As we move through the super- 
lattice from the VO2 layers to the TiC>2 barrier layer 
we see the DOS evolve from the typical bulk BVO layer 
form (layers 1 to 3) to a new form near the heterojunc- 
tion. In particular, we see a strong peak in DOS in layer 
4 near the Fermi energy Ep with 1.6 e~ per V ion in 
agreement with our conjectures. As far as magnetic or- 
der is concerned (Fig. [9]) , we find G-type AFM order in 
which neighboring transition metals ions have opposite 
spins both within and between layers for the bulk like 
VO2 layers (layers 1 and 2). The magnetic orders in 
layer 3 and 5 are reconstructed to be ferromagnetic with 
the magnetic moment 0.9755 (ib per V ion and the mo- 
ment 0.4750 fiB per Ti respectively. The magnetic order 
in layer 4 is ferrimagnetic with a net magnetic moment 
0.2170 hb- In the corresponding band structure illus- 
trated in Fig. [ITjJ we identify three partially filled two- 
dimensional bands that cross the Fermi level, and sixteen 
full d-bands. These bands accommodate 9 d-electrons 
obtained by taking the average d-charges per transition 
metal ion and summing the charges of all the layers, as 
illustrated in Fig. [T2] discussed below. The two majority 
spin partially-filled bands, and the single minority-spin 



FIG. 8: (color online) Layer-resolved d-projected density of 
states at Ti and V sites versus energy. 



partially-filled band, are responsible for the extra charge 
at V sites in layer 4 and 6 as shown in Fig. [9l The three 
full minority-spin bands, which are lower in energy than 
the partially-filled minority-spin band but closer to the 
Fermi level than other full minority-spin bands mainly 
contribute to minority-spin charge concentrated in lay- 
ers 4, 5, and 6. The wave functions for those bands 
are strongly hybridized within those layers. The lowest 
majority-spin band in Fig. [10] is mainly of V character 
in layer 4 but lower in energy due to the more attractive 
local potential at this site. The peak of the DOS for the 
minority spin at the V site seems to shift higher toward 
the Fermi energy than expected. This is due to the fact 
those states are strongly hybridized with the states at 
the Ti sites which are shifted toward higher energy, an 
effect not included in the qualitative analysis. (See the 
potential profile of the realistic conjecture of case (b) in 
Fig. QJ. Other full bands are responsible for BVO bulk 
states away from the region in which the electronic inter- 
face reconstruction occurs. 

In Fig. QTJ we plot plane-averaged Hartree potentials. 
The first layer starts with VO2 at z = 0. Transition 
metal layers repeat with a spatial period of 3.96 A with 
the center (z = 15.84 A) layer occupied by Ti02 layer. 
At z = 1.98 A, the first BaO layer starts and repeats 
also at a period of 3.96 A; the local minima closest to 
Ti02 correspond instead to LaO layers. By inspecting 
the local minima in the Hartree potential, we can see 
that the DFT calculations largely verify the scenario of 
case (b) in Fig. [T] except that the dashed line wiggles in 
the region of bulk BVO. The wiggling may due to the 
slight difference of the strength of ionic bonding between 
Ba-0 and V-0 along the layer direction, which causes the 
small dipole modification of charge distribution in BVO 
ionized layers. 

In Fig. 1121 we show the corresponding electron lin- 
ear number density A(z) = 0.5 J p(x,y, z)dxdy along the 
layer-growth direction. This charge number spatial dis- 
tribution p{x, y, z) include the total charges contributed 
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FIG. 9: (color online) Spin and layer resolved densities-of- 
states versus energy. The labels on top are the layer labels 
indicated in Fig. [2] Since we allow G-type AFM order in the 
barrier material, there are two distinct transition-metal ions 
in each layer. The DOS for majority spins and minority spins 
are represented by red and blue lines respectively. 




Momentum 

FIG. 10: (color online) Spin-resolved band structure. The 
horizontal axis represents superlattice crystal momentum. 
The red solid lines represent the majority-spin bands while 
the blue lines represent minority-spin bands. The momentum 
dependence along the layer direction fc z is captured in the re- 
gion between the F and Z points. The three partially occupied 
bands associated with EIR are flat along k z , demonstrating 
their quasi-two-dimensional character. 



FIG. 11: (color online) Plane averaged Hartree potential 
along the layer-growth direction z. The solid dashed line indi- 
cates the net potential variation between charge-concentrated 
centers at each layer. 
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FIG. 12: (color online) Linear electron number density X(z) 
along layer-growth direction z. The dashed lines represent the 
boundaries between different layers used in constructing the 
layer decomposition of this charge density. The numerical 
numbers between dashed lines in each layers are the total 
number of electrons from d-derivative bands in Fig. 1101 In 
this case, the total number of electrons is nine. 



from the ideal conjecture we made for this case and al- 
most the realistic conjecture. In general, the charges in 
the well are transferred from the TiC>2 layer to adjacent 
VO2 layers as discussed in section II. 



from the bands between and - 0.9 eV with respect to 
Fermi level in Fig. [TDJ We see an apparent charge ac- 
cumulation (1.5144 e~) at the TiC>2 layers centered at 
z = 11.88 A and z = 19.8 A along with much less charge 
(0.6295 e~) at the TiC>2 layer. Remote transition metal 
layers have small deviations from a d 1 valence. The total 
amount of charges from the (i-bands is 9 e~ by adding 
the charges at each layer in agreement with our predic- 
tion in case (b). The calculation shows small deviations 



V. SUMMARY AND DISCUSSION 

In this paper we have explored some possibilities for 
the realization of artificial two-dimensional electron sys- 
tems by growing quantum wells consisting of a few atomic 
layers of one layered oxide surrounded by another ma- 
terial. In particular we consider systems which involve 
heterojunctions between polar insulating perovskites like 
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LTO and nonpolar insulating perovskites like BVO. Per- 
ovskites can be viewed as consisting of alternating AO 
and BO2 layers where B is a transition metal. The main 
idea of this paper is that inversion symmetric quantum 
wells will necessarily have a different number of AO and 
BO2 layers, and that this will necessarily lead to an 
electronic reconstruction which induces two-dimensional 
electron systems near the quantum well when one ma- 
terial is polar and one is nonpolar. Four different cases 
can be identified, in which the quantum well material is 
either polar or non-polar and in which it either has an 
excess AO layer or an excess BO2 layer. The excess lay- 
ers induce cither a unit increase or a unit decrease in the 
total (Z-orbital occupancy per transition metal site when 
summed over all layers. This charge must be localized 
close to the quantum well to avoid large band shifts due 
to polarization driven electrostatic potentials. As long 
as inversion symmetry around the quantum well center 
is not-broken spontaneously, half of this charge must ap- 
pear on each half of the quantum well center. Because 
of the strong-correlation character of the oxide electronic 
structure, these charges will tend to be localized in half- 
filled LSDA+U bands which are weakly hybridized so 
that both cross the Fermi energy. 

We have tested and elaborated on this qualitative pic- 
ture of the electronic structure of polar-nonpolar hetero- 
junction quantum wells in perovskites by performing self- 
consistent LSDA+U electronic structure calculations for 
a representative system composed of BVO and LTO with- 
out any atomic relaxation at the interface. The LSDA+U 
calculations capture the interplay between long range po- 
lar electrostatic energy and local physics driven by ener- 
getic offsets between different transition metal ions and 
by p — d hybridization. We find that in the two cases we 
have examined one additional two-dimensional bands is 
induced by the energetic offsets. In all cases the quan- 
tum well structure creates a two-dimensional metal inside 
a three dimensional Mott insulator, and alters the char- 
acter of the magnetic order. The reconstructed magnetic 
order is ferromagnetic in the active well region and an- 
tiferromagnctic in bulk region. The appearance of two- 
dimensional fcrromagnctism in intimate exchange con- 
tact with antifcrromagnetism is intriguing from the point 
of view of spintronics. 

Because of the complexity of oxide materials, our cal- 
culation should be viewed as an exploration of possi- 
bilities rather than as predictive in any detailed sense. 
Among the many sources of uncertainty, the Hubbard U 



interaction used in our LSDA+U is a phenomenological 
parameter for which we have chosen a standard bulk val- 
ues. This value should likely be renormalized near the 
quantum well. More challenging to account for is the 
role of inevitable lattice distortions, oxygen stoichiome- 
try variations, and other defects near the interfaces. The 
interfaces between the thin film oxide and the substrate 
on which it is grown can also play a complicating role 
if either the barrier material or the substrate is polar. 
For a few cases we have explored the influence of lattice 
distortion effects by allowing the atoms to relax along 
the layer growth direction in the active region, maintain- 
ing Uv=Uri= 3.2 eV as before. What do we find is that 
the amount of charge transfer to or from particular layers 
and the details of the magnetic orders are highly sensitive 
to ionic relaxation. These complications make the accu- 
rate prediction of the electronic reconstruction and the 
associated magnetic order highly nontrivial. It is also de- 
sirable to understand if exotic symmetry breaking phases 
such as superconductivity may appear in the interface . 
Further progress in understanding detailed properties of 
the two-dimensional electron systems we envision which 
almost certainly require inputs from experiments. The 
microscopic calculations reported on here can be used as 
guidance in constructing phenomenological models which 
can be fit to any data which might emerge from future 
work and serve as the theoretical basis for possible exotic 
phases in the interface. The present paper has addressed 
d 1 transition metal systems. Considerations similar to 
those explained here also apply to transition-metal oxides 
with heavy transition metal elements (Cu, Ni) in which 
the t2 g bands are fully occupied far below the Fermi level, 
and low-energy excitations are expressed within the e g 
bands. 
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